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We derive the expressions for the 1-loop corrections in cosmological, Eulerian, perturbation theory 
to the matter bispectrum and to the qalaxy bispectrum, assuming local galaxy bias, in presence of 
non-Gaussian initial conditions. We compute them explicitly for the particular case of non-vanishing 
initial bispectrum and trispectrum in the local model and for a non- vanishing initial bispectrum alone 
for the equilateral model of primordial non-Gaussianity. While the primordial contribution to the 
matter bispectrum for values compatible with CMB observations is dominant over the component 
due to gravitational instability at large scales, 1-loop perturbative corrections due to non-Gaussian 
initial conditions correspond to just a few percent of the gravity-induced bispectrum at mildly non- 
linear scales, similarly to what happens for the matter power spectrum. However, in the perturbative 
expansion for the galaxy bispectrum, 1-loop diagrams arising from non-linear bias are responsible for 
significant large-scale contributions, indeed exceeding the primordial component, both for the local 
and equilateral model. We study the peculiar dependence on scale and on the shape of the triangular 
configurations of such additional terms, similar in their origin to the large-scale corrections to the 
halo and galaxy power spectra that raised significant interest in the recent literature. 



I. INTRODUCTION 

Recent observations of the Cosmic Microwave Background (CMB) have confirmed that primordial perturbations 
in the matter density are well described by a Gaussian random field, whose statistical properties are encoded in its 
two-point correlation function. Constraints on higher-order correlations such as the three-point function or its Fourier 
transform, the bispectrum, are model-dependent but generically consistent with a level of non-Gaussianity below 
0.1%, [TJ[5]. Indeed, the simplest models of inflation typically predict departures from Gaussianity a few orders of 
magnitude smaller then current limits, Oil]. It is therefore evident that a possible detection in future observations of 
a significantly larger non-Gaussian component would rule-out a single-field slow-roll model for inflation and greatly 
improve our ability to discriminate among the currently wide range of available phenomenological models of the early 
Universe (for a review see, for instance, [S]). 

Measurements of the bispectrum of the CMB temperature fluctuations have been, so far, the most powerful probe of 
primordial non-Gaussianity. On the other hand, the effects of non-Gaussian initial conditions on the growth of matter 
perturbations and on large-scale structure observables have been studied for quite a long time, focusing, in particular, 
on corrections to the high-mass tail of the cluster abundance [SJ |71 [HI ISl HHJ Ell dl 1131 [HI [E] and to cumulants 
[TBI [T7\ [T8] and higher-order correlation functions [19, 20, 2TJ [22] of the galaxy distribution. Until recently, the 
galaxy bispectrum has been considered the most promising observable, expected to provide comparable or even better 
constraints on non-Gaussianity from measurements in future high-redshift, large- volume galaxy surveys, [211 123] . 

This picture radically changed during the last year after Dalai et al. [23 showed, in numerical simulations, that 
a relatively small departure from Gaussianity can have a large effect on the power spectrum of dark matter halos. 
These results, somehow anticipated by early theoretical works on the correlators of the distribution of peaks in the 
matter density field ^ , attracted significant attention [271[2HJ[21[3ni[3Il[Ml[Ml[31[MllMll371ISHllMl[ini[lI], 
finding further confirmation in independent numerical results, |30L 1341 155] . Most importantly, they allowed to derive 
constraints on the amplitude of primordial non-Gaussianity of the order of CMB limits from current large-scale 
structure observations, [271155] . 

These first works focused their attention on the effects on the power spectrum of biased objects, halos and galaxies, 
proposing quite different theoretical descriptions. For instance, starting from the expression defining the specific non- 
Gaussian model assumed for the N-body simulations of Dalai et al. [24], a correction to the halo bias can be derived 
in the framework of the peak-background split, [211 [271 [33] . This approach, however, seems to be hardly extensible to 
different models of primordial non-Gaussianity. An alternative description, proposed by Matarrese and Verde |32| and 
based on the statistics of high-peaks in the matter density, provides a correction to the power spectrum of the peak 
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distribution in term of the generic primordial matter bispectrum and can be therefore applied to any non-Gaussian 
model. A similar result, but derived in terms of 1-loop perturbative corrections to the galaxy power spectrum under 
the assumption of local, Eulerian bias in position space, has been suggested by Taruya et al. [37]. Finally an even 
different interpretation has been considered by McDonald [33^. 

While all the proposed corrections to the linear bias function of the halo power spectrum present the same functional 
dependence on redshift and scale, which is quite generic when non-Gaussian initial conditions are considered, they 
also make distinct and relevant approximations. Despite the fact that a satisfactory agreement with the halo power 
spectrum measured in numerical simulations has been shown by Grossi et al. ^30 when a correction factor, possibly 
accounting for ellipsoidal collapse of virialized objects, is included, we believe that more work is still required to 
provide a proper theoretical model for this effect applicable on a large range both in mass and scales and to a generic 
model of primordial of non-Gaussianity. 

A natural test of the different approaches mentioned above consists, in this respect, in the derivation of predictions 
for higher-order correlation functions, and, in the first place, for the bispectrum. A first step in this direction has 
been made by Jeong and Komatsu [31], which applied to the bispectrum the theoretical description of Matarrese and 
Verde '32] based on the distribution of high density peaks. In this paper we pursue the same goal extending instead 
the results of Taruya et al. [37\ for the matter and galaxy power spectra to the matter and galaxy bispectra. For the 
galaxy bispectrum, we therefore assume a local bias prescription and compute all tree-level and 1-loop contributions 
in Eulerian Perturbation Theory up to a given order in the linear matter density field. We will show that some of these 
contributions, including a term corresponding to the one considered by Jeong and Komatsu [3_lj, are indeed dominant 
at large scales, in a similar way as in the power spectrum case. We will show, moreover, that, unlike the power 
spectrum, the galaxy bispectrum receives significant corrections also for an equilateral model of non-Gaussianity. 

It is reasonable to assume, based on our results, that previous predictions for the constraints on primordial non- 
Gaussianity from measurements of the galaxy bispectrum in current and future redshift surveys, [201 121] are likely 
to improve significantly. In fact, we can conservatively expect that for a non-Gaussian parameter, as for other 
cosmological parameters, the bispectrum should be able to provide at least comparable results as the ones provided 
by power spectrum measurements, W3. A proper assessment of the expected constraints on non-Gaussianity, 
however, will have to wait future work, and, in particular, necessary comparisons with N-body simulations. It is worth 
reminding that our assumption of a local bias does not include a prediction for the corresponding bias parameters, 
which are usually derived, for Gaussian initial conditions, in the framework of the Halo Model. However, it should 
provide the functional form of the most relevant contributions to the galaxy bispectrum which is the main result of 
this paper. 

In Section [H] we discuss relevant inflationary models responsible for a possibly observable non-Gaussian component 
in the density perturbations. We introduce in particular the local and equilateral expressions for the bispectrum of the 
initial curvature perturbations that will be used as examples in the rest of the paper. In Section [ITT] we summarize the 
Perturbation Theory (PT) approach to gravitational clustering of the matter density field and review previous results 
on the matter power spectrum in presence of non-Gaussian initial conditions. We then derive the 1-loop corrections 
to the matter bispectrum and the tree-level contributions to the matter trispectrum, i.e. the four-point function 
in Fourier space. In Section |IV| we consider the perturbative expressions for the correlators of the galaxy density 
field under the assumption of local bias. In particular we derive the 1-loop corrections to the galaxy bispectrum and 
compute them for the local and equilateral models of primordial non-Gaussianity. In Section [V| we discuss similarities 
and differences of our results with those recently proposed by Jeong and Komatsu [31J. Finally, we present our 
conclusions in Section [VT] 



Initial conditions for the evolution of matter perturbations are usually given in terms of the early-times correlators 
of the matter density field, 5, or of the curvature fluctuations, (f>. For Gaussian initial conditions, all statistical 
properties of the initial fields are encoded in their two-point function, or the power spectrum, in Fourier space. Small 
departures from Gaussianity can be characterized by non-vanishing three- and four-point functions, or their Fourier 
transform, respectively, the bispectrum and trispectrum. In the following sections we will study the effects on the 
late-times matter and galaxy bispectrum due to two common phenomcnological models for non-Gaussian, primordial 
curvature perturbations, the ZocaZ model [22] [44] [45| [46] [47] , and the equilateral model, [48] . 

The local model describes a small correction to Gaussian initial conditions, represented by the expression for 
Bardeen's curvature perturbations <i>, in position space 



II. MODELS OF PRIMORDIAL NON-GAUSSIANITY 
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where (j) represents the Gaussian component and the second and third term on the r.h.s. the quadratic and cubic, non- 
Gaussian components. We assume here the CMB convention for the definition of the f^L parameter, which implies 
that the curvature $ is evaluated during matter domination. The local model typically describes non-linearities in the 
relation between the inflaton and curvature perturbations |44U46llTf] . and in general models where non-Gaussianity is 
produced outside the horizon such as curvaton models |49j , inhomogeneous reheating [50. .51j or multiple field inflation 
[52]. 

From Eq. ([T]) one can derive the following expression for the leading contribution to the curvature bispectrum, 

Si°=(A:i,fc2,fc3) = 2/ArLF$(A:i)F$(fc2) + 2pcrm., (2) 

with the curvature power spectrum P$(fc) defined in terms of the Gaussian component alone as ((/)(ki)(/)(k2)) = 
(ki2)P$(fci), where we introduce the notation = k^-l-kj. Similarly, the curvature trispectrum will be given by, 

4-(ki,k2,k3,k4) = 4/2,^P$(fci)P$(fc2)[P*(fci3) + i^$(fci4)]+5perm. 

+ Q9NLP^{ki)P^{k2)P^{k^) + 3 perm. (3) 

Notice that this quantity depends on six variables, which can be chosen to be the magnitudes of the four wavcnumbers 
plus the magnitudes of two sums like, for instance, ki2 and ki^. 

The factor /^^ in front of the first term on the l.h.s. of Eq. ([s]) is clearly a consequence of the definition of Eq. (jlj. 
However, inflationary models predicting an initial bispectrum corresponding to Eq. (|2|, can generically predict an 
initial trispectrum of the form of Eq. ([s]) but with a different normalization factor, sometimes denoted as tat^ (see, 
for instance, |53j). Such distinct quantity typically presents a dependence on the parameters of the inflationary model 
different than Jnl- While we will consider, for simplicity, Eq. (|3| as our model for the initial local trispectrum, we 
should keep in mind that the relation between the amplitudes of 5$ and T$ is, in general, less trivial. 

Current limits on the value of from measurements of the CMB bispectrum correspond to —4 < /j^^- < 80 
(95% C.L.), [1], while large-scale structure constraints, exploiting the effect on the power spectrum of highly biased 
objects, are given by —31 < < 70 (95% C.L.), [36], although somehow different analyses provide higher values, 

m- 

As a complementary expression to the bispectrum of the local model, Babich et al. |48] introduced an equilateral 
model of primordial non-Gaussianity defined by the bispectrum 

+Pl'\ki)Pl'\k2)P^{k^) + 5 perm.] . (4) 

This functional form assumes larger values for equilateral triangular configurations, that is when ki k2 — k^, while 
the local model is more sensitive to squeezed configurations where fci ~ fc2 <C ^3. The normalization of the expression 
in Eq. Q has been defined in a such a way as to provide an identical value as the local expression, Eq. ^ for 
equilateral configurations, ki = k2 = k^ when = f^j^. It has been shown, [351|S31|5S], that the equilateral model 
for the bispectrum closely approximates the bispectrum predicted by inflationary models with non-canonical kinetic 
terms like DBI inflation [51?, Ghost inflation [57 or higher derivatives [58l[59]. Current limits on the amplitude of 
equilateral non-Gaussianity are —151 < < 253 (95% C.L.), [2 . 

While we limit ourselves to these two models of primordial non-Gaussianity, it should be stressed that several other 
functional forms for the initial bispectrum with peculiar dependencies on the shape of the triangular configurations 
are predicted by infiationary models in the literature (see [55 and references therein). 

The matter overdensity in Fourier space (5k is related to the curvature perturbations <i>k by Poisson's equation as 

Skiz) = M{k, z) $k, (5) 

where we introduce the function 

with T(fc) being the matter transfer function and D{z) the growth factor. Notice that since we conform to the 
CMB convention for the definition of the non-Gaussian parameter /tvl, we denote with <1> the 'primordial curvature 
perturbations, i.e. evaluated during the matter dominated era, not their value linearly extrapolated at present time. 

The initial correlators of the matter density field are related to the correlators of the curvature perturbations, in 
general as 

(<5ki • • • <5k„ ) = Af (fci , z) • • • Af (fc„ , z) ($ki • • • *k J , (7) 



4 



so that, in particular, the Hnear power spectrum is given by 

Poik)^ M^k,z)P^{k), (8) 
while the initial bispectrum and trispectrum are given respectively by 

Bo{ki,k2,k3) = M{ki)M{k2)Mik3)B^{ki,k2,k3), (9) 

and 

To(ki,k2,k3,k4) = M(A:i)M(fc2)M(fc3)M(A:4)r*(ki,k2,k3,k4). (10) 

In the following sections we will study the effects of the local and equilateral models of primordial non-Gaussianity 
on the evolved matter bispectrum and trispectrum in the framework of Perturbation Theory and, with a similar 
perturbative approach based on the assumption of local galaxy bias, on the galaxy bispectrum. 

III. MATTER CORRELATORS AND NON-GAUSSIAN INITIAL CONDITIONS 

In Fourier space, the perturbative solution to the equations of gravitational instability for the matter overdensity 
can be formally expressed as [60] . 

S.^si'^ + S^^^+Sl'^ + ..., (11) 



where 



(") 



where we keep implicit the time dependence as S^'' ~ D(z), D{z) being the linear growth factor, and where 
_F'„(qi, ...,q„) is the symmetrized kernel for the n-th order solution. 



As shown explicitly in [61], the expansion of Eq. (Ill leads to an ill-defined perturbative solution for the power 
spectrum, characterized by large cancellations between contributions of the same order. New approaches to the 
physics of gravitational instability have been introduced with the Renormalized Perturbation Theory in fM] 1621 163] , 
and with the Renormalization Group formalism in |64j leading to a well behaved perturbative expansion of the matter 
power spectrum and more accurate predictions of simulations results. A first step in extending these techniques 
to higher-order correlators can be found in [BS]. We will leave the extension of renormalized perturbation theory 
or renormalization group approaches to non-Gaussian initial conditions to future work (but see Izumi and Soda 
[66]), limiting ourselves to observe that traditional perturbation theory can still provide a satisfactory description of 
the impact of non-Gaussianities on the matter bispectrum, particularly since, as we will show, this turns-out to be 
negligible with respect to the Gaussian component. We also notice that, with this work, we are filling a gap in the 
literature that could serve as a starting point for future investigations. 

For all numerical evaluations in this Section, as well as in the following, we assume a flat ACDM cosmology with 
the WMAP 5-yr, [5], cosmological parameter values = 0.044, = 0.214, ffo = 71.9, = 0.96, Ut = 0.0, r = 0.087 
and with a r.m.s. of the matter fluctuations corresponding to CTg — 0.8. 



A. The matter power spectrum 



In this section we summarize, for completeness and because they will be useful in the following sections, previous 
results in cosmological perturbation theory relative to the matter power spectrum. Specifically, 1-loop perturbative 
corrections to the matter power spectrum have been studied in [67' ^681 for Gaussian initial conditions. Additional 
contributions up to the same order in PT for generic non-Gaussian initial conditions have been derived in [37. and 
computed for the local and equilateral non-Gaussian models assuming a non-vanishing initial bispectrum. Notice 
that, while in the Gaussian case all 1-loop corrections are of the same order, O(<5o), in terms of the linear matter 
overdensity Jg, for non-Gaussian initial conditions, we have an extra term — denoted as P12 in [37 — of order 0((5q). 
Similarly, 2-loop corrections consist of contributions of several different powers in the initial field 6q. We choose to 
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write the perturbative solution for the matter power spectrum as a series of powers of the initial field Sq, indicating 
each contribution of order 0{5q) as so that 

P{k) = P^^^k) + P(3)(fc) + p(4)(^) + 0(^5)^ (^3) 

where, for Gaussian initial conditions, all terms with odd index vanish. Introducing the notation 

6D(k^2)P^,ik^) ^ (<5«4f ) + eye, (14) 
we have P^^^ = Pn = Pq being the linear matter power spectrum, while 

p(3)(fc) = Pi2, (15) 

pW(fc) = Pi, + pU + Pi, + P{i, (16) 



with 



Pi2 = 2 yd3gP2(q,k-q) Po(fc,g,|k-q|), (17) 

P4 = 2 Jd^qF^{q,k-q) Po{q) Po(|k - q|), (18) 

Pii = y"Arf'pJ^2(p,k-p)P2(q,-k-q) To(p,k-p,q,-k-q), (19) 

P/3 = 6 Po(fc) yAP3(k,q,-q) Po(g), (20) 

2 y"d^gdV3(P,q,k-p-q) ro(-k,p,q,k-p-q). (21) 



pll 
^13 



In the equations above Bq and Tq represent the initial bispectrum and trispectrum, respectively, related to the higher 
order correlators of the curvature perturbations according to Eq. ([9]) and Eq. (10). P12 is the 1-loop non-Gaussian 



contributions studied by [37], the terms Pfj and P/3 constitute the Gaussian 1-loop corrections, |S71inH], while Pj, 
and P/3^ are 2-loop corrections vanishing for Gaussian initial conditions. We will consistently ignore, for simplicity, 
in our evaluations throughout the paper all 2-loop contributions. 

The contributions to the solution for matter correlators in cosmological perturbation theory can find a representation 
in terms of diagrams (see f5D]), often helpful to identify different terms and their characteristics. In Fig. [l] we show the 
diagrams representing the tree-level, 1-loop and 2-loop contributions to the matter power spectrum in the expressions 
above. The upper half shows, in particular, the elements entering the diagrammatic representation, given by verteces 



corresponding to the perturbative solution of Eq. (Ill and the initial correlators joining them. A simple line connecting 
two points corresponds to the linear power spectrum Po{k) while the three connected lines — as those in the P12 term — 
correspond to the initial bispectrum Bq and the four connected lines — as in P22^ and P(^ — correspond to the initial 
trispectrum Tq. It is easy to see how non- vanishing higher-order correlators such as Bq and Tq provide several extra 
contributions to the non-linear matter power spectrum with respect to those present for Gaussian initial conditions, 
corresponding to a larger variety of possible diagrams. This will be even more evident for the non-linear bispectrum. 



B. The matter bispectrum 



We consider now the non-linear matter bispectrum for non-Gaussian initial conditions. Similarly to power spectrum, 
we write the perturbative solution for the matter bispectrum up to corrections 0{Sq) as 

P(fcl, fc2, h) ^ P(3) ^ ^(4) ^ ^(5) ^ 5(6) ^ 0(^7)^ (22) 

where P^^^ ~ 0('^o) where all odd contributions vanish for Gaussian initial conditions. Here P^'^-' — Pm = Po 
corresponds to the initial bispectrum, while higher-order terms are given by the following individual contributions 

p(4)(fci,fc2,A:3) = Pfi2 + 5n2, (23) 
p(^)(fci,fc2,A:3) - 5(22 + 5(22 + S(i3 + Bn3, (24) 
p('^)(fci,fc2,A:3) = P^22 + 5(23 + 5(23 + 5n4 + 2-loop terms, (25) 



This is different from tfic more common notation wfiere p(p) indicates the p-loop correction. 



Diagram elements: verteces and initial correlators 



F (r, n n \ 





Matter power spectrum: tree— level, 1—loop and 2— loop diagrams 



^12= • 




+ cyc. P^22 = 




^13 ~ 





+ eye. 



FIG. 1: Elements of the diagrammatic representation of perturbativo corrections of matter correlators and diagrams representing 
the tree-level {Pii), 1-loop (P12, P22 £md Pl-^) and 2-loop {P22 and P13) contributions to the matter power spectrum up to 
terms of order 0{6o). Lines connecting two points correspond to the linear power spectrum Po{k), the three connected lines in 
the P12 term corresponds to the initial bispectrum Bo and the four connected lines in P22^ and PH correspond to the initial 
trispectrum To. 



where we make use of the notation for the individual terms on the r.h.s. of the equations above 

SD0^123)BiJl{h,h,ks) = (4? 4?) + eye, 

anafogous to the one for the power spectrum. We can distinguish the tree-level expression 

Bfi2 = 2 i^2(ki,k2) Po(fci) Po{k2) + 2 perm., 

from the 1-loop corrections 

Bn2 = y^rf^9 (q, ks - q) To (ki , k2 , q, ka - q) , 



5(22 = 2Po(A:i) 



P2(ki,k3) J^2(q,k3-q) Bo(fc3,g,|k3-q|) + (fc3^fc2) 



2 perm. 



= F2(ki,k2) [Po(fci) Pi2{k2) + Po{k2) Pi2(fci)] + 2 perm., 



(26) 



(27) 



(28) 



(29) 
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Bii2 = 4yd3qF2(q,k2-q) F2(ki+q,k2-q) So(fci,g,|ki+q|) Po(|k2-q|) + 2pcrni., (30) 

Bfi3 = Wo{ki,k2,k3) y"AF3(k3,q,-q)Po(g) + 2perm., (31) 

Bii^ = 3Po(fci) Jd^q F3(ki,q,k2 - q)Po(A:2,<?, Iks - q|) + {ki ^ ^2) + 2 perm., (32) 

^222 = 8 ydV2(-q,q + ki)P2(-q-ki,q-k2)F2(k2-q,q)Po(g)Po(|ki + q|)Po(|k2-q|), (33) 

= 6Po(fci) y"d3gP3(ki,k2-q,q) P2(k2-q,q) Po(|k2-q|) Po(g) + 5perm., (34) 

B'lk - 6 Po(fci) Po(A:2) P2(ki,k2) J d\ P3(ki,q,-q) Po(9) +5 perm. 

= P2(ki,k2)[Po(A:i) Pi3(A:2) + Po(fc2) Pi3(fci)]+2perm., (35) 
Pfi4 = 12 Po(fci) Po(fc2) jd\ P4(q,-q,-ki,-k2) Po((z) + 2 perm.. (36) 

Again, we ignored 2-loop corrections to the P^^-* term. Here, the 1-loop contributions PI22! ^123; ^123 ^^"^ ^114 
are those derived by Scoccimarro [69' for Gaussian initial conditions, while the others, P(x2' ^122' ^1221 -^113 

depend on the initial bispectrum Pq and trispectrum Tq. Initial correlation functions of order higher than the 
trispectrum appear only in 2-loop corrections and higher. In Fig. |2]we show the diagrams representing the tree- level 
and 1-loop contributions to the matter bispectrum described by the expressions above. 

In Fig. [3] we show the equilateral configurations of the matter power spectrum P(fc, k,k) as a function of k for 
non-Gaussian initial conditions. We choose values for the non-Gaussian parameters within current 1-cr uncertainties. 
In the case of the local model we include the effect due to the initial trispectrum to 1-loop corrections, while for the 
equilateral model only a non-vanishing initial bispectrum is considered. The upper panels show the effect of local 
non-Gaussianity with = 50 and qml = 10"* at redshifts z = 1 {left) and z = 3 {right). The fourth-order term 
P^''^ represented by long-dashed line, includes the correction due to the initial trispectrum Pffj of Eq. (28 1 for the 
local case, which however, is completely negligible, so that the long-dashed line effectively indicates the tree-level term 
B[i2 alone. The unnaturally large value for the parameter g^L is assumed only for illustrative purposes, since only 
for such large values we obtain a contribution comparable to the one from the term proportional to in Eq. 
In the lower panels of Fig. [3] we consider the equilateral model with /]^^ = 100 at redshifts z ~ 1 {left) and z = 3 
{right). 

The P*^-^) = Pq contribution, represented by short-dashed lines, is, by definition of the equilateral model itself, 
identical in the local and equilateral model for equilateral configurations, however the P*-^) contribution is in principle 
different as it corresponds to integrals over a continuum of different configurations of the initial bispectrum. As it 
is the case for the power spectrum studied in [3^, the extra 1-loop contributions corresponding to P'-^^' are quite 
suppressed with respect to the Gaussian component. Specifically, in the range of scales where the P^^-' term is 
largest, roughly corresponding to 0.03 to 0.05/iMpc~^ in our cosmology, it represents a correction of a few percent 
to the Gaussian prediction, still smaller then the primordial contribution P^'^^ by nearly one order of magnitude for 
equilateral configurations. Notice that, for the equilateral model such configurations present the largest value for 
the primordial component, while this is not the case for the local model. Finally we notice that Gaussian 1-loop 
corrections in the P*-^-* term are significant, as expected, at somehow smaller scales, above 0.1/iMpc~^ at redshift 
z = 1 (see |7ni |7T] for recent simulations results with Gaussian initial conditions) . 

To study the dependence on the triangle shape, it is convenient to introduce a reduced matter bispectrum [72] . 
defined as 

r^fh h 1 \ - B{ki,k2,k3) 

Q{ki,k2, fcs = , 37 
S(fci, fc2, fca) 

with 

S(/ci, fc2, ^3) = P(fci) P(fc2) + 2 perm. (38) 



Its perturbative expansion in powers of 60 can be derived from Eq. ( 22 1 and from the series 



S(fci, fc2, ^3) - + + + 0{5l), (39) 



8 



Matter bispectrum: tree-level diagrams 




Matter bispectrum: 1-loop diagrams 




FIG. 2: Diagrams representing the tree-level and 1-loop contributions to the matter bispectrum with generic non-Gaussian 
initial conditions. 



with 

S(4) = p(2)(^^)p(2)(^^)^2perm., (40) 

perm., (41) 

j;(6) _ p(3)(/j^)p(3)(fc2) -I- 2 perm. -h p(2)(A;i)p(^)(fc2) + 5perm., (42) 

so that, with the same convention for the notation, we obtain 

Q{kiM, h) = Q^-'^ + + Q^'^ + Q^'^ + 0(^o)- (43) 
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FIG. 3: Equilateral configurations of the matter bispectrum, 
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fSh = 100 at redshifts z = \ (left) and z = 3 (right). The mail 
the primordial component B^^K 
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B{k,k,k), for non-Gaussian initial conditions. Upper panels: 
(left) and z = 3 (right). Lower panels: equilateral model with 
I effect due to non-Gaussian initial conditions is represented by 



where 



r(3) 

<?'-" ^ fjj, (44) 



B(4) _ Q(-1)I](5) 

B(6) „ Q(l)l](5) „ g(0)s(6) _ Q(-1)E(6) 



Q^°^ = (45) 
Q^'^ - ^""""w.r"' ^ ^ (46) 

^ (47) 



The Q*-^^' and Q*-^' contributions vanish for Gaussian initial conditions. 

In Fig. |4] we show the reduced matter bispectrum Q{ki,k2,k3) with fixed ki and k2 = 1.5 fci as a function 
of the angle 9 between ki and k2 for non-Gaussian initial conditions of the local kind with /nl = 50 and non- 
vanishing initial bispectrum and trispectrum^ . The upper panels assume ki = 0.01 /iMpc~^, while lower panels assume 
ki = 0.05 ft,Mpc~^. For the first set of configurations the leading correction to the Q'^^^ term {long-dashed line) is given 
by the primordial component Q^^^-* (short-dashed line), which assumes the largest values for 9 ~ n, corresponding 



^ Again, the contribution depending on the initial trispectrum To, even assuming gpfL = 10*, is completely negligible for all configurations 
shown in Fig. [4| 
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FIG. 4: Reduced matter bispectrum Q{ki, k2, fcs) with fixed ki and fc2 = 1.5fci as a function of the angle 9 between ki and k2 
for local non-Gaussian initial conditions with /jvl — 50 and qnl = 10*. Upper panels assume fci = 0.01 /i Mpc~^ while lower 
panels assume 0.05 /iMpc~^. Predictions in the left and right panels are evaluated respectively at redshift z = 1 and 

2 = 3. As in Fig.lsl Q^~^^ (short-dashed line), which includes the primordial component B^^\ represents the main non-Gaussian 
contribution, which assumes the largest values for 9 ~ n, corresponding to the squeezed limit for this set of configurations. 



to the squeezed limit for this set of configurations. For the second set of configurations with ki = 0.05/iMpc~^ the 
leading correction to Q^^^ is given instead by the 1-Ioop contributions in the term Q^^^^ (dot-dashed line), larger at small 
values of 9 corresponding to large values of k^. Predictions in the left and right panels are evaluated respectively at 
redshift z = 1 and z = 3. Choosing a value for the ratio ^2/^1 closer to 1 would lead to a much greater non-Gaussian 
component in the squeezed limit since, in fact, Q^~^\k,k,k^) ~ l/fc| and therefore diverges for 6 ^ it. We avoid in 
our examples such extreme cases of squeezed triangles since configurations with fci = ^2 3> ^3 are difficult to measure 
in actual surveys due to the limited resolution in the wavenumber value. 

Fig. |5] shows the same results as Fig. |4] but for the equilateral model with /^vl = 100. In this case the largest 
non-Gaussian corrections correspond to the nearly equilateral configurations with 9 ~ O.Ttt. 

From these two examples it is particularly evident the different effect that alternative models of primordial non- 
Gaussianity can have on the matter bispectrum and that could lead, in the optimistic case of the detection of a large 
signal in future observation, to distinguish the specific shape-dependence of the initial component and to possibly 
identify the corresponding model. 



C. The matter trispectrum 



We consider here, finally, the expression for the matter trispectrum up to the same order in perturbation theory 
considered above for the bispectrum, that is 0(<5q), since it will be useful in the next section. We will focus our 
attention to the tree-level contributions, neglecting loop corrections. Following the notation introduced for the power 
spectrum, we can write 



T(ki, k2, k3, k4) = t(4) -I- t(5) + r(6) + &(5l), 



(48) 
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and, denoting individual terms as 

<5D(kl234)r.,7™(kl,k2,k3,k4) ^ (<S^^ ^'.^ ^3 ) + CyC, (49) 

we have, in the first place, the initial trispectrum T^^^ = Tim = To, i.e. the connected four-point function, vanishing 
for Gaussian initial conditions, then 

T(5)(ki,k2,k3,k4) = Tim + Tlli^ (50) 

with TIii2 being the tree-level correction due to gravitational instability and a non-vanishing initial bispectrum Bq 
given by 

Tlm{kuk2, fcs, k^) = 2 Bo(fci, fc2, A^ia) [i^2(ki2, k3)Po(A:3) + F2(ki2, k4)Po(fc4)] + 5 perm., (51) 

while T(Ii2, which we do not write down explicitly, represents a 1-loop correction involving the connected five-point 
function, and finally 

r(<')(ki,k2,k3,k4) = r/i22 + Tim + 1-loop and 2-loop terms, (52) 

where 

r/i22(fci,fc2,fc3,fc4) = 4^^2(ki3,-ki) F2(ki3,k2)Po(fci)i"o(A:2)Po(fci3) + llperm., (53) 
Tlmiki,k2,k3,h) = 6 P3(ki, ka, k3)Fo(fci)Po(fc2)-Po(fc3) + 3 perm., (54) 

represent the leading tree-level contribution for Gaussian initial conditions. Again we neglect all 1-loop and 2-loop 
corrections depending, in general, on the square of Bq, on products of Pq and Tq and on higher-order initial correlation 
functions. 
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Matter trispectrum: tree-level diagrams 




+ 11 perm. 




+ 11 perm. 



FIG. 6: Diagrams representing the tree-level contributions to the matter trispectrum with generic non-Gaussian initial condi- 
tions. 



In Fig.|6]we show the diagrams representing the tree-level contributions Tim, 7iii2, 2Y122 ^^^d T/^^g to the matter 
trispectrum as described above. 

We will show in the next section that the primordial component To = Tini and the Tl^^2 term, arising from a 
combination of primordial non-Gaussianity and non-linearities in the matter evolution, will contribute significantly 
to the 1-loop corrections to the galaxy bispectrum. As we do not consider, for the time being in this paper, all 2-loop 
corrections we ignore consistently 1-loop corrections to the matter trispectrum as they will be part of contributions 
to the galaxy bispectrum involving one further integration, therefore effectively corresponding to 2-loop corrections 
in terms of the linear density field. 



IV. GALAXY CORRELATORS IN THE LOCAL BIAS MODEL AND NON-GAUSSIAN INITIAL 

CONDITIONS 



Taruya et al. |i37j derived an perturbative expression for the galaxy power spectrum in presence of non-Gaussian 
initial conditions following the local bias prescription of Fry and Gaztafiaga [73]. They have shown that an effect 
on the linear halo bias similar to the one measured in N-body simulations fM] [SHI IMl ISH] in its redshift and scale 
dependence can be derived as a 1-loop correction to the galaxy power spectrum which is of order 0{5q) in the linear 
density field 5q. Here we follow the same approach, based on the assumption that relation between the smoothed 
galaxy overdensity (5g(x) and the smoothed density contrast (5ij(x) is local, that is 

5gA^)^f[5R{^)]. (55) 

Assuming small fluctuations at large scales one can Taylor-expand the relation above as 

^,.«(x) = b,6n{^) + y [4(x) - (4(x))] + ^[4(x) - (4(x)] + (56) 
where hi are constant bias parameters. In Fourier space, for A: ^ 0, we can write the filtered galaxy overdensity as 

OO J 

= E (57) 



where 



4?A(k) = d\i ■ ■■d\n5D{\^-q.i...n)5R{cii) ■ • •'5fl(qn), (58) 
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with 5_R(k) = Wii{k)6]i, Wr being an appropriate smoothing function. From a comparison between these expressions 
and Eq.s ( |TT| ) and ([T2| , it is evident that perturbative corrections to the correlators of the galaxy distributions can 
be obtained in a similar way as those for matter correlators considered in the previous section. 

In this section, we compute the galaxy bispectrum up order 0((5q) in perturbation theory, excluding therefore 
1-loop corrections to the bispectrum for Gaussian initial conditions — which are of order 0{Sq) — but including 1- 
loop corrections involving the initial bispectrum, and, for the local non-Gaussian model, the initial trispectrum. This 
amounts to neglect non-linear effects at small scales while concentrating our attention on the large-scale contributions. 
For the filter Wr we assume a top-hat window function in position space with R = 5/i~^Mpc, although this value, 
used for the quantities shown in the figures, is not strictly relevant to the main results of the paper. 



A. The galaxy power spectrum 

We begin by summarizing the results of Taruya et al. |37j on the galaxy power spectrum in the framework of the 



local bias model of Eq. ( 56 ) and non-Gaussian initial conditions that will be useful for our analysis of the galaxy 
bispectrum. Following Smith et al. |74j and Taruya et al. |37j . we define the unfiltered galaxy power spectrum Pg{k) 
in terms of the power spectrum of the smoothed density field Pg^nik) simply as Pg{k) = P^{k) /W'^{kR). The 



perturbative expression for the galaxy power spectrum, consistent with the expansion of Eq. (56 1, can be written as 



Pgik) ^ ^^^^ = PgM + Pa,12 + Pg,22 + Pg,13 + ■■■, (59) 



where each term on the r.h.s. is defined according to the notation 

Sn{k,,)Pg,,{k,) = wik,R)W{k,R) ■ ^'^^> 

We choose therefore to denote each contribution to the galaxy power spectrum, and later to the galaxy bispectrum, 
in terms of their dependence on the bias parameters bi, instead of their perturbative order, keeping in mind, however, 
that, in general, Pg^ij contains terms 0((5q'^"') and higher. The first contributions, involving corrections up to 0((5g), 
can be written in terms of the non-Zmear matter correlators P, B and T as |37l 1741 175] . 

Pg^n = blP{k), (61) 
Pg,i2 - 61^2 /a W^fl(q,k-q) S(fc,g,|k-q|), (62) 



Pa.22 = I Id'q M^^(q,k-q) P(g) P(|k-q|), (63) 



^/.22 = ^ Jd\d^pWR{ti,k-ti)WR{p,^k-p)T{q,k-q,p,-k-p), (64) 
PgS3 = bib3 P{k) Jd^q W^q) P{q), (65) 
Pgis - ^ ydV'pl^fl(q,p) W^fl(q+P,-k-q-p) r(k,q,p,-k-q-p), (66) 



where we introduced, for convenience, the combination of window functions 



M/ t \ W{q,R)W{q2R) 

l^^(qi,q2)^ W^(,,,P) • 



Clearly, the quantities Pg.22 and Pg^is as defined in Eq. (60 1 correspond to the sums Pg 22 + -^"922 and Pg 13 + Pg\d.^ 
respectively. It is as well evident that the diagrams in Pig. [T] can represent as well the terms in the perturbative 
expansion of the galaxy power spectrum once we replace the mode couplings in each node with the powers of the 



series in Eq. (56 1 and the primordial correlators with those of the non-linear matter density field. 

We can rewrite the expressions above making explicit the perturbative contributions to the matter power spectrum, 
bispectrum and trispectrum to keep track of the order of the expansion in terms of the initial density field 5q. Excluding 
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corrections of order 0{Sq) we have 



Pg.ll 



P' 



pll 



-^3,13 



pll 



p(2) _^ p(3) ^ p(4) 



&162 Jd^qWR{q,k-q) [B^^\k, q,\k-q\) + B^^\k, q,\k~q\)], 
I Jd^qd^pWR{q,k-q) W^fl(p,-k-p) (q, k-q, p, -k-p), 

jd\d^pWR{q,p) W?«(q+p,-k-q-p) rW(k, q, p, -k-q-p). 



Ms 
3 



(68) 
(69) 

(70) 
(71) 
(72) 
(73) 



For consistency with the previous section and for simplicity we will ignore all 2-loop contributions. The interesting 
term in this expansion is given by the 0(<5q) correction, 

P%{k) = 6162 (d\ WR{q,k-q) Bo{k,q, \k-q\), 



(74) 

which depends on the initial matter bispectrum B^"^^ = Bq. In fact, assuming the expression for the initial bispectrum 
of the local model given by Eq. ^ , in the large-scale limit A: ^ we have 

P^g%{k) = 2fNLbib2 ld\WR{q,k-q) M{k,z)M{q,z)M{\k-q\,z) 
X [P^{k)P^iq) + P$(fc)P<,(|k-q|) + F$(g)P$(|k-q|)] 



'^^^ AfNLbib2M{k,z)P^ik) jd\Wl{q) M\q,z)P.^{q) 
= 4.fNLbib2M{k,z)P^ik)aUz), 



which can be written as 



corresponding to Eq. (5.3) in |37j . Excluding corrections of the order 0((5^), the galaxy bispectrum is given by 



pAk)^p^'h+p^%+p^-:L 



(3) 



5(3) 



tr>\ 

and, neglecting the small term Pg = bfPi2^ we can write 



Pg{k) C hi 



M{k,z) 



Poik), 



(75) 



(76) 



(77) 



(78) 



and the Pg,i2 contribution can be expressed as a scale-dependent, non-Gaussian correction to the linear bias, defined 
as 



blifNL) = b.ifNL = 0) + Ab,ifNL,k), 



and given by 



A61 = 2/]^£- 62 a^(2) 



M(fc, z) 



(79) 



(80) 



As pointed-out by Taruya et al. [37 , this expression is similar in its scale and redshift dependence to the correction 
derived in the peak-background split framework by |24l |27l 13^ , from a model for high-peak correlators by [32 and in 
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an even different derivation, but related to the local bias assumption, in [35]. In particular, the peak-background split 
approach, [24l|27l|36] provides, for the correction to the linear bias of halos of mass M, bh,i{M, z), the expression 

Ah,,, = 2fSl[hhAM,z) - ^]5cj^y (81) 

where 5c — 1.686 is the critical threshold for spherical collapse. 

Following a different approach, Matarrese and Verde [32] considered the peak number density defined by the 
expression pp^Af(x) — 9[5r{^) — 5^, where is a step function. Based on earlier works by Grinstein and Wise [25] 
and Matarrese et al. [26] which provide an expression for the peak two-point function of the form 

Cp(xi,X2)=exp[X(xi,X2)]-l, (82) 

where the X(xi,X2) represents a series involving all higher-order correlators of the initial density field, they derive 
a correction to the peak correlation function in terms of the three-point correlation function. Their result can be 
rewritten in the form of a peak linear bias correction given by, in our notation 

A^i = ^ 1^^^ J d\ WR{ci, k - q)i3o(fc, q, |k - q|) (83) 

where v = 5c/ (Jr- This expression is valid in the high-peak limit v/ur 3> 1 but, as the one derived by Taruya et al. 
|37j . is quite general since it applies to any form of the initial bispectrum Bq. In the large-scale limit and for local 
non-Gaussianity we can find the asymptotic behavior for the, Lagrangian, peak bias correction 



where the last equality shows the equivalence with Eq. (81 1 after identifying hp^i — 1 + v/ctr as the Eulerian linear 
peak bias. 



It is interesting to notice that, if we interpret the bias parameters introduced by the local bias model of Eq. (56 1 
as halo bias functions, bh.i{M), and if we assume for these quantities the Gaussian predictions of the Halo Model, 
[76|, 177] then, in the high-threshold limit v/ctr S> 1, we can approximate 



so that all these different bias corrections, including Eq. (80 1, provide similar quantitative results. In fact, the factor 
^/(t|j in Eq. (]83|, could just as well be interpreted as a non-linear, quadratic bias parameter. 



Taruya et al7^37 notice how the presence of the cr|j term in the asymptotic expression of Eq. ( 80 1 might represent 
a problem due to its logarithmic diverge as i? — > in the case of local non-Gaussianity. However, if we are to assume 
that indeed such different results describe the same effect, the peak-background split approach would provide a well 
defined prediction for the product b2<j\ when 62 describes quadratic non-linearities of the halo bias relation. On the 



other hand, Eq. (81 1 has been shown to well describe numerical results, [30] |34l [38] , up to a corrective factor q = 0.75, 
possibly related to ellipsoidal collapse [3D]. They also notice that other difficulties arise for negative values of 



or 62, since in those cases, in the large-scale limit, the galaxy power spectrum in Eq. (78) becomes negative. Such 
problem, however might just be an artifact of the truncated perturbative expansion rather than a pathology of the 
local model of non-Gaussianity. 

We will not discuss further all peculiar approximations involved in the different approaches proposed in the literature. 
As we will se in the next Section, large-scale behaviors proportional the quantity ur are present for several corrections 
to the galaxy bispectrum. We will therefore avoid interpreting our results as quantitative predictions limiting ourselves 
to present them as a derivation of the functional forms that we can expect for non-Gaussian contributions to the galaxy 
bispectrum. 



B. The galaxy bispectrum 

We study now the perturbative corrections to the galaxy bispectrum under the assumption of local bias, which 
represent the main results of this paper. Similarly to the case of the power spectrum, we can define the unfiltered 
galaxy bispectrum as 



(86) 
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which we can expand as 

^9 = -^3,111 + -^3,112 + BgS22 + -Bg,113 + • ■ • , (87) 

where each term can be formally defined as 

^ (\c \B (h h h\-^^ ^' (Ji,>i)'^j;(k2)'^yk3)) + perm. 

The contributions containing terms up to 0{Sq), as a function of the non-linear matter correlators P, B and T are 

Bg^ni = h\B, (89) 
Blii2 = blb2 W^fl(ki,k2) P(fci) P(fc2) +2 perm., (90) 

Sg"ii2 = ^ yA W?ii(q,k3-q) T(ki,k2,q,k3-q) + 2perm., (91) 
Sg^22 = ^ W^i?(ki,k3) P(fci) M?ji(q,k3-q) P(fc3,'7,|k3-q|) + 5 perm., (92) 

Bls22 = bibl jd\ W?fl(q,k2-q) W?i^(ki+q,k2-q) B(fci, g, |ki + q|) P(|k2-q|) + 2 perm., (93) 

T 2 L /• 

^3^13 = ^ W^fl(ki,k2) P(fci) yA M?fl(q,k2-q) P(fc2,g,|k2-q|) + 5 perm., (94) 

Sg"ii3 = ^ &?&3 B{k^M.h) jd\ WUq) P{q). (95) 

^3^,222 = f yAW^fl(-q,q + ki)W^fl(-q-ki,q-k2)M?fl(k2-q,q)P(g)P(|ki+q|)P(|k2-q|), (96) 

Bg,i23 = ^y^^yi^(klk2)P(A;l) jd\wU^M~ ci)P{\k2~ ci\)P{q) + b x^eru,., (97) 

^9^123 = ^^i&263Wjj(ki,k2)P(fci)P(fc2) J qW^q) P (q) + 2 peim. , (98) 

<ii4 = ^W^fl,(ki,k2)P(A;i)P(fc2) jd\wl{q)P{q) + 2^evT^., (99) 

where we ignored 2-loop contributions in the terms i?g,i22 and Bg n^ which depend on the matter five-point function 
and 2- and 3-loop contributions in the terms Pg.222, Bg.123 and Bg m depending on squares of B, products of P and 
T and on the matter six-point function. 

The 1-loop corrections in Bg 2221 Bg 123' -^9^123 and Bg ^ only depend on non-Gaussian initial conditions via small 
corrections to the matter power spectrum at mildly non- linear scales, [37. , and have therefore no effect on the large-scale 
behavior. Following McDonald [75 and Smith ct al. [74 , the terms Pg'i23 and Bg -^^^^ can be interpreted as providing a 
redefinition — or renormalization — of the bias parameters since they present the same functional dependence as Bg ^^2- 
Also, the not-renormalizable terms Pg 222 and Bg 123 behave, for equilateral configurations, like P^{k) as fc — > 0, and 
we can expect them to be sub-leading with respect to a significant primordial contribution to the galaxy bispectrum. 

Since we are mainly interested in the effects due to primordial non-Gaussianity at large scales, we will compute 
all perturbative corrections to the galaxy bispectrum up to order &{Sq) in the linear matter density field. So- We 
will therefore ignore the 1-loop corrections to the matter bispectrum corresponding to the P'^-* term as well as the 
non-linear bias terms corresponding to Pg.222, ^9,123 and Bg 114. Introducing the notation 

Bg^ijk{ki,k2,k3) = bibjhkBg,ijk{ki,k2,k^), (100) 



which simply factorizes the bias parameters from the quantities defined in Eqs. (89 1 to (99 1, we can rewrite the 
remaining contributions in terms of the expansions of Eq.s (13), (22 1 and (48 1 for the matter power spectrum P, 
bispectrum B and trispectrum T, respectively, to keep track of the perturbative order in terms of the linear matter 
density Jq. We obtain, for the Pg.iii contribution simply the expansion of the matter bispectrum of Eq. (22 1, that is 

Pg 111 = P(3) + ^(4) + b(5) + 0(^6)^ (^0^) 
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FIG. 7: Contributions to the galaxy bispectrum with local non-Gaussian initial conditions assuming /jvl = 50 and gjvi ~ 10* 
and unitary bias parameters, fei = 62 = 63 = 1 at redshift z = 1: i3c/,ii2 {upper left panel), Bg,i22 {upper right panel) and -85,113 
{bottom left panel). See the text for a detailed explanation of all terms. The bottom right panels shows a comparison of all 
different contributions, Bg.in {continuous black line), 59,112 {dashed blue line), Bg,i22 {dotted magenta line) and Bg.ns {dot- 
dashed green line) with the thin lines corresponding to the values of Bg,iii and i3g,ii2 for Gaussian initial conditions. -85,122 
and -Bg,ii3, excluding corrections 0(5o) are vanishing in the Gaussian case. Assumes a smoothing scale oi R = 5h~^ Mpc. 



where each component has been shown in the upper panels of Fig. [3] For the -Bg, 112 contribution we have instead two 
terms of order 0{Sq) and two of order 0{Sq), so that we can write 

Ba,ii2^Bl%, + B%, + 0{S^,) (102) 

where 3^%, = B^j^ + B^^iit^ with 

Bg,n2 = WR{ki,k2) P^^Hki)P^^Hk2) + 2peYm., (103) 
b"{[^2 = ^y^'g W?i?(q,k3-q) T(4Hki,k2,q,k3-q) + 2perm., (104) 

and where B^^^ = B^'iH + bI^{[1^ with 

Bifil = W^i?(ki,k2) p(2)(fci)p(3Hfc2) + 5perm., (105) 
Bl!n2 = ^yd=^gH?fl(q,k3-q) T(5)(ki,k2,q,k3-q) + 2perm.. (106) 

All this terms are show, for non-Gaussian initial conditions of the local kind and assuming f^L = 50 and qnl = 10"*, 
in the upper left panel of Fig. It] Specifically, -Bg 112 {short-dashed line) represent the usual, Gaussian tree-level 
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prediction for galaxy bispectrum induced by quadratic non-linearities in the bias relation, while i?g'|x2 {dotted line) 
is a small correction to the same quantity due to the 1-loop contribution to the matter power spectrum P^^^ for 
non-Gaussian initial conditions. 

More interesting are the terms derived from Bg\i2 depending on the matter trispectrum. For the local model we 
consider the non-vanishing initial trispectrum of Eq. ^ given by the sum of a term proportional to and one 
proportional to gNL- We plot these two components separately as ^g^ii2''/jv^ (medium- dashed line) and Bg\[^^g^^ 
( long-dashed line) , with Bg^{[^ = Bg^-[[^j:^^-\-Bg^{[^g^^ . One can see that these terms provide the largest contribution 
to Bg, at large scales, with a particularly strong scale dependence. The large value gNL = 10** as been chosen here 
for illustration purposes only, in order to provide a contribution comparable in size to the one proportional to f^L- 
In other words, for more natural values of gpf^ ~ few, we can expect this term to be negligible. Assuming the 
term to be the leading contribution to the initial curvature trispectrum given by Eq. ^ in the local non-Gaussian 
model, we have 

^"n2!/„J^i>fc2,fc3) ^ 2f^^ Jd'qWR{q,k,-q) M{k,)M{k2)M{k:,-q)M{q) 

X {P*(fci)P<t,(fc2) [F*(ki3-q) + P<i,(ki + q)] + 5 perm.} + 2 perm. (107) 
For equilateral configurations, in the large-scale limit this expression can be approximated as 

^"n2l/„.(fc.fc>fc) M^ik) Pl{k) jd\ Wl{q) M\q) P^{q) 

= 36 fljL al M\k) Piik). (108) 

This asymptotic behavior can be compared to the one of the primordial contribution, B^"^^ = Bq, given, for equilateral 
configurations, by Bo{k,k,k) — 6 fNLM^{k)P^{k), so that 

B"/Z.Ak,k,k) ,_.o 4(z) D{z) 

Boik,k,k) - ^^^^M(M ^ ^ 

where we made explicit the dependence on redshift. We can therefore expect this component to be dominant at 
large scales, but its relevance with respect to the primordial one, decreasing with redshift. Notice that the redshift 
dependence is in part determined by the uj^iz) factor which might find a physical interpretation, as discussed at the 
end of the previous Section, when considered in its product with the quadratic bias parameter 62- For our choice of 
the smoothing scale, R = 5h^^ Mpc, we have (Jr{0) = 1.07 and we keep this factor in all evaluations shown in the 
figures. 

We can also compare the asymptotic behavior with the Gaussian, tree-level prediction for the matter bispectrum, 
B1^2j which is of the same order in PT. We have, again for equilateral configurations, 

B"{[tf„Ak,k,k) ,^o7 aliz) ^ 

Bi,2ik,k,k) ~ 3 ■'^^M^k,z)^ k^' ^ ' 

where the ratio shows an even larger scale dependence, but is independent of redshift, as expected. 

It is interesting to notice that the fifth-order contribution B^^^^^ depending on the fifth-order matter trispectrum 



correction T^^^^ of Eq. (106 1, vanishing for Gaussian initial conditions is also equally important. We show this 



contribution as the dot-dashed line in the upper-left panel of Fig. [7] where one can see how it provides a significant 
contribution at intermediate scales between the trispectrum induced, 1-loop correction B^^^y^ and the Gaussian tree- 

this term arise due to the coupling between the initial bispectrum 



level term Pg_'ii2- As discussed in Section 



IIIC 



and quadratic non-linearities in PT. From the expression for Tl^^yi Eq. (51 1 one can write 

K'M^kxMM) = Jd^qWR{q,k3-q) 

X {Bo{ki, fc2, fci2) [P2(ki2, k3-q)Po(k3-q) + F2(ki2, q)Po(g)] + 5 perm.} 

4^2 perm., (Ill) 

and in the large-scale limit, but for generic configurations and generic initial conditions one can show that 

^i'm (^1,^2,^3) yPo(fci,fe,fc3)a|+ ^Po(fti) Jd\WUq)Bo{k2,q,q)+perm., (112) 
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FIG. 8: Same as Fig.[7]but for equilateral non-Gaussianity with fNL ~ 100. In this case only a non- vanishing initial bispectrum 
is considered. 



which shows how one of the permutation in Eq. (51) provides a term proportional to initial matter bispectrum Bq, 
while the others an integration of the same. This means that we can expect this contribution to retain in part the 
shape-dependence of the initial bispectrum. Assuming local non-Gaussianity, for equilateral configurations, we can 
compare Bg\[2^ to Bq, obtaining the expected scale-independent behavior at large-scales 

K!n2ik, k,k) fc^o 262 
Bo{k,k,k) " ^ ' 

which shows as well how the relative weight of this contribution decreases with redshift. The ratio with the Gaussian 
term B{i2 is given instead by 

Bi'Al\k,k,k) fe_.o 131 aliz) D{z) 

Bi^2{k,k,k) 3 ■'^^M{k,z) fc2 ■ ^ ^ 

In Fig.|8]we show the same results as Fig.[7]but for the equilateral model, assuming /^vl = 100- Notice that in this 
case we consider only a non- vanishing initial bispectrum Bq, with no contributions from the initial trispectrum Tg 
of the matter field, although they should be in principle taken into account for any realistic model. For this reason 
the Bg^{['^ terms are not present. It is remarkable, on the other hand, that the Bg^{[^ contribution due to the 
term of Eq. (51 1 in the expression for the matter trispectrum — dot-dashed line in the upper left panel of Fig. |8] — still 
provides a large correction at large-scales. This behavior is specific to the galaxy bispectrum as it has been shown 
that lowest order correction to the galaxy power spectrum for acceptable levels of equilateral non-Gaussianity are 
ncghgible and probably unobservable, [29l |37j . 

We are left to consider the Bg 122 and Sg,ii3 terms contributing to our series for the galaxy bispectrum up to 
0(^q). These are shown in the upper right and bottom left panels, respectively, of Fig.jTjfor local non-Gaussian initial 
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conditions with fj\jL = 50. Excluding corrections 0((5g), both -63,122 and i3g, 113 are vanishing for Gaussian initial 
conditions. 

The Bg 122 consist in turn of the two terms, both of the fifth order in 

K% = ^ W^K(ki,k3) p('Hfci) W?/?(q,k3-q) -B(3)(fc3,q,|k3-q|) + 5perm., (115) 

^g"i22^ = yAH^fl(q,k2-q) M^fl(ki + q,k2-q)B(^)(fci,-?,|ki+q|) P(2)(|k2_q|)+2perm.. (116) 

Bg 122 is represented in the upper right panel of Fig. [rjfor the local model and of Fig. [sjfor the equilateral model by 
the dotted magenta line, while Bg\22 by the dashed line. The continuous line is their sum. Their asymptotic behavior 
at large-scales, for equilateral configurations and assuming local initial conditions is given by 

Bl%{k, k, k) 12 fNLP^ik, z)^^^ - D{zf [loc. NG], (117) 

Bl'i^{k, k, k) 12 fNLP^iK z)^j^Jd\ W^{q)P^{q) ^ [loc. NG], (118) 



M{k,z) 

where we remark a peculiar form of scale-dependence in Bg^{^2 i while for equilateral initial conditions, their contri- 
bution is suppressed at large scales 

Blfl{k,k,k) 12 fNLM^{k,z)Pl,'\k) jd\ WUq)M\q,z)Pl'\q) ^ D{zf k^ [eq. NG], (119) 

Bl'Al\k,k,k) \2 fNLM{k,z)Pl'\k)jd\wi{q)M\q,z)Pl'\q)P^{q)^ D{zf k [eq. NG]. (120) 

Finally, the Bg contribution is made up by the two terms 113 and 113 shown, respectively by the dotted 

and dashed lines in the bottom left panel of FigJTlfor the local model and of Fig. [8] for the equilateral model. It is 

I I J h ^ / — ^ I 

clear from a comparison of Eq. (94 1 with Eq. (^92| that term P^ 113 is identical to -Sg'122 discussed above and it is 

therefore relevant at large-scales only for local non-Gaussianity. The term B^^^fr^ instead is simply contributing to 
the overall multiplicative factor of the initial bispectrum, since 

^gfif3^ = ^5o(fci,fc2,A:3)4W. (121) 



The term Bg of Eq. (94) therefore provides contributions equal to the -Bg.iii term to all orders in the perturbative 
expansion of the matter bispectrum. However, in our choice to retain only corrections up to 0(i5q), the sole contribution 
proportional to the primordial bispectrum B^^^ = Bq is present. This is consistent, in a way, with focusing our 
attention to the large-scale contributions due to primordial non-Gaussianity. 

We remind that, in this approximation, the only relevant initial correlators are the bispectrum Bq and the trispec- 
trum Tq since we ignored as well the contribution of the initial five-point function to the trispectrum term T*^^\ 
in principle consistent with this expansion, but leading to an effective 2-loop correction with one integration in the 
expression for T^^-* itself and a second one in the expression for Bg\i2 above. Taking into account the identical 
contributions — up to a multiplicative factor — we could therefore rewrite Eq. ( 87 1 as follows 



b\ [l + ^^4) 4.111 + b\b2 (bI,,2 + Bl'112) + ^? (I + I) 4' 122 + ^^1^25^122, (122) 



which shows that, in our approximation, only five different functional forms are present. Each multiplicative factor 
in the expression above should be interpreted in terms of the proper — renormalized — bias parameters to be compared 
with measurements in numerical simulations, |74l 1751 . 

To explore the shape-dependence of the galaxy bispectrum, as in the case of the matter bispectrum, we can define 
a reduced galaxy bispectrum as 

n (u I. u \ - Bg{ki,k2,k3) 

Qg[kl,k2, fc3 = ^ , , ^ , (123 

Eg(fci,fc2,fc3) 

where 

Sg(fci,fc2,fc3) = ^9(^1)^9(^2)+ eye. (124) 
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Similarly to the bispectrum we introduce the notations 

^g,ij,kii^i^^'2^^3) = Pg,ij{ki)PgM{k2) + perTon-; (125) 

= bibjbibk^g^ijM{ki,k2,k3), (126) 

Qg,2]k{ki,k2,k3) = QgAjk{ki,k2,k3), (127) 

to keep track of the non-linear order in the bias expansion and factorize the bias parameters. We can consider then 
the approximation 

Qg{kl, ^2, fca) — J^Qg.lll + ^Qs,112 + p^Qg,122 + ^Qg,113 (128) 

where we exclude higher-order contributions corresponding to corrections of order 0((5g) and where 

g, (129) 



^9,111 



^9,112 



(130) 



BgS12 — Qg,lll'^g,ll,12 



Bg.l22 — Qg, 112^3, 11, 12 " Qg.lll ( ^9,11,22 + Sg,i2,12 
^9,122 = — ^ ^ (131) 

A -B9413 - Qg,lllS^g,lia3 .-,oo\ 
^9,113 - • (,13^j 

^9,11,11 

In Fig. [9] we plot the different contributions to the reduced galaxy bispectrum with local non-Gaussian initial 
conditions, assuming &i = 62 = 63 = 1, as a function of the angle between the vectors ki and k2 with fc2 = 1.5 fci for 
two different values of the magnitude ki and at redshifts z = 1 {left panels) and z — i {right panels). Fig. 10 shows 
the same results for equilateral non-Gaussian initial conditions. Thick lines represent the sum of all contributions to 
each term Qg^ijk while the thin lines for Qg^m and (5g,ii2 correspond to their values for Gaussian initial conditions. 

In these plots it is particularly evident the difference between the two non-Gaussian models, despite the fact that 
we did not choose a set of triangular configurations including particularly squeezed triangles since, at most, for = tt 
we have ^2 = 1.5 A:i = 3 k^. Still, one can observe for the local model in Fig. [9] a significant increase in the non- 
Gaussian signal as — > TT, both in the correction to the matter bispectrum contained in Qg.iii and even more so in the 
lowest order correction due to non-linear, quadratic bias Qg,ii2. Indeed, as expected from the plots of the equilateral 
configurations of Fig. [71 Qg,ii2 is responsible for the largest non-Gaussian signature. Clearly, as mentioned before, 
choosing a set of triangles with ki ~ k2 would lead to an even greater effect, and, in the limiting case ki — k2, the two 
contributions diverge as — *■ tt. Such extreme triangles are, however, difficult to measure in real surveys and we limit 
ourselves to provide an idea of the non-Gaussian corrections to the galaxy bispectrum for quite generic triangular 
configurations. The signature of an equilateral primordial non-Gaussianity is, on the other hand, evident, in the plots 



of Fig. 10 where the largest contribution is expected for 9 ^ 0.7 tt, the closest to an equilateral configuration. Again, 
regardless of the effect due to the initial trispectrum which we do not include here, the largest effect is on the Qg,ii2 
term, which departs significantly from the constant predicted by the Gaussian, tree-level expression, due essentially 



to the -Bgfii2'' contribution discussed above, that is 



QgM2 = ^ for Jnl = (133) 

which corresponds to the value 1 in the plots where we set &i = 62 = 63 = 1. Indeed we find that such correction, for 
a level of non-Gaussianity well within current constraints such as /nl — 50, can be a factor of a few larger than the 
Gaussian prediction already at fc ~ 0.01 /iMpc~^ and at redshift z = 1, that is, at scales nearly accessible by current 
surveys. The overall impact on the galaxy bispectrum depends, of course, on the relative size of non-linearities in 
the galaxy bias relation, or, in simple terms, on the ratio &2/&1- Such ratio, assuming the halo model predictions for 
Gaussian initial conditions, is expected to be small, and in fact negative, for galaxy samples like the main samples of 
the Sloan Digital Sky Surveys (SDSS) or of the 2dF Galaxy Redshift Survey, but it is also expected to increase with 
an increasing value of the linear bias, b\, and be therefore particular significant for populations like the Luminous 
Red Galaxies sample of the SDSS, {TSU^HD]- 
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FIG. 9: Galcixy bispectrum with local non-Gaussian initial conditions, assuming /jvl ~ 50 and gjvL = 10 , as a function of the 
angle between the wavenumber ki and k2 with k2 — 1.5 ki for the different values 0.01 {left panels) and 0.02 /i Mpc~^ {right 
panels) and at redshift z — 1 {upper panels) and z = 3 {lower panels). Assumes unitary bias parameters bi — b2 = bs — 1 and 
a smoothing scale oi R = 5 Mpc. 



We can also notice that, while the non-Gaussian effect on the mostly due to the primordial component 

Bq to the matter bispectrum, increases with redshift, the effects due to the combination of non-linear bias with 
non-Gaussian initial conditions instead decrease with redshift. This observation is in contrast with the notion that 
high-redshift surveys will provide the best constraints on non-Gaussianity, based in turn on analyses which took into 
account only the detectability of the primordial component, |21j . Quantitive predictions for the expected constraints 
on non-Gaussian parameters from measurements of the galaxy bispectrum in upcoming surveys will however require 
further work and we leave them entirely for future publications. 



V. ON THE HIGH-PEAKS THREE-POINT FUNCTION 



In this Section we comment on the differences between our results and the similar ones derived by Jeong and 
Komatsu [31 , following a different approach that extends the study of the high-peak two-point correlation function 
and power spectrum of Matarrese and Verde |32l to the high-peak bispectrum. 



As mention already at the end of Section IV A Matarrese and Verde fT2^ derived a correction to the peak two-point 
correlation function ^p,M in terms of the three-point correlation function. 



Cp(xi,X2) ~ — Co,i?.(xi,X2) + — Co,i?.(xi,X2,X2), 



(134) 



'R 



where ^o,fl f^nd C,q^r correspond to the two- and three-point function of the initial density perturbations, filtered on a 
scale i? and where v = 5c/<jr with being the r.m.s. of the smoothed linear fluctuations. 
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From our assumption of local bias, Eq. (56 1, we can write a similar expression for the galaxy two-point function, 

^g(xi,X2) ~ 6iC_r(xi,X2) + 6ife2Ci?,(xi,X2,X2), (135) 

where and (^r are the filtered and e^o/iierf matter two- and three-point functions evaluated at the redshift of interest. 
Both expressions give rise to similar corrections to the power spectrum, particularly large when the initial three-point 
function assumes large values for squeezed triangular configurations, as it is the case for local non-Gaussianity. We can 
notice, in fact, that the three-point functions that appear in Eq. (134) and (1351 are evaluated for collapsed triangles. 



The second expression, however, includes all higher-order corrections to the matter bispectrum in PT as described by 
Eq. (22 1. In addition, a comparison between the two leads to the improper identifications bi ^ v/or and 62 ~ I'^'r 
for the bias parameters in the high-peaks limit, 3> 1. In fact, these values correspond to the high-threshold limit 
for the Eulerian expressions of the halo bias parameters derived in the framework of the ellipsoidal collapse, [771 181j . 



but represent Lagrangian quantities in Eq. (1341. 

Matarrese et al. [26], on which the results of Matarrese and Verde 32J are based, derived as well an expression for 
higher-order correlators and, in particular, for the three-point function we have, 

Cp(xi,X2,X3) = F(X1,X2,X3) [^p(xi,X2)Cp(xi,X3),^p(x2,X3) 

+fp(xi,X2)^p(x2,X3) + 2 perm. + ^p(xi,X2) + 2 perm. + 1] 
-^p(xi,X2) + 2 perm. - 1, 



(136) 



where 



{00 n-2n-j-l I I \n 

n=3j = l k = l ■'' ■' '' 



(137) 



with ^p"^ representing the smoothed linear matter n-point function with the n arguments given by j times 

xi, fc times X2 and {n — j — k) times X3. Expanding the exponential for small values of the correlators to include 
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corrections up to the matter four-point function, one obtains 

"3 3 ly^ 



F(X1,X2,X3) ~ 1 



3 Co,k(xi,X2,X3) + ^^^|!j^]j(xi,X2,X3,X3). 



(138) 



'R 



Retaining this order of corrections, proportional to the forth power in the Hnear density field, Jeong and Komatsu 
[31j derive from Eq. (136 1, the expression for the bispectrum of high-peaks given by 

,,3 



Cp(xi,X2,X3) 



1 ,(4) 

'R 



3-Co.J?(xi,X2,X3) + - —Co A(''i'^2,X3,X3) + 2 perm. 

R Z CTe ' 



-Co,fl(xi,X3)^o,K(x2,X3) -f 2 perm. 



(139) 



Again, we can compare this expression with the similar one that can be obtained from the local bias prescription 
of Eq. (|56|) , that is 



{Sg{^i)Sg{^2)Sg{^3)) - 6?(^(xi)5(x2)5(x3)) + hlb2{S{^i)S{^2)S\^3)) + 2 perm 



(140) 



where the four-point function ((5(xi)5(x2)(5^(x3)) is given by a connected and a non-connected component correspond- 
ing, respectively, to the second and third terms on the r.h.s. of Eq. (1391, so that the galaxy three-point function is 
given by 



Cff(xi,X2,X3) 



&iCfl(xi,X2,X3) -I- -6i62Cfl^(xi,X2,X3,X3) -I- 2 perm. 



+fe?&2Cfl(xi,X3)^_R(x2,X3) +2 perm. 



(141) 



In the approach of Section IV however, the correlators (ji a nd (or their unsmoothed yersions) are fully evolved, 
while the corresponding quantities ^o,-Ri Co,R. Siiid ^g^jj of Eq. (1391 are their initial counterparts. This implies that 
Eq. (139) does not include non-Gaussianities due to gravitational instability but only those determined by non- linear 
bias and the initial conditions. Jeong and Komatsu |31j . however, assume that the iniiiaZ correlators correspond to 
relatively late-time {z ~ 10) quantitities so that, in their expression for the halo bispectrum, based on Eq. (1391, 
the matter correlators do include contributions due to the growth of structures. Their result is therefore formally 
equivalent to the one of Eq. ( |141 1, leading to the same non-Gaussian corrections, particularly since, under their 
assumption, the predictions of Eq. (139) for the bias factors loose their validity and are replaced by parameters to be 
fitted by comparison with observations. 



VI. CONCLUSIONS 



In this work we studied several aspects of the effects of non-Gaussian initial conditions, in the framework of Eulerian 
perturbation theory, both on the matter as on the galaxy bispectrum under the assumption of local bias. 

In Section[TTT] we have computed perturbative 1-loop corrections to the matter bispectrum in presence of a primordial 
non-Gaussian component corresponding to two common phenomenological models resulting in a non- vanishing initial 
matter bispectrum, and, for the local model of non-Gaussianity, also in an initial matter trispectrum. We have shown 
that, as it is the case for the power spectrum, the effects of such corrections, given current constraints from CMB 
observations on the amplitude of non-Gaussianities both of the local and equilateral kind, at mildly non-linear scales 
are of the order of a few percent of the leading contribution to the matter bispectrum due to gravitational instability. 

In Section |IV| we considered then a perturbative expansion from a local bias prescription and the corresponding 
corrections to the galaxy bispectrum, and identified several large contributions enhanced by non-linear bias that 
might exceed the direct effect due to the initial component. In particular we have shown that, for local primordial 
non-Gaussianity, a significant non-Gaussian effect is expected due to the non-vanishing matter trispectrum, both by its 
primordial component in the local model, depending on the initial trispectrum, at very large scales, thus confirming 
recent findings by Jeong and Komatsu |31j . as well as by the next-to- leading order correction in perturbation theory, 
depending instead on the initial bispectrum, which appears to be larger at smaller scales. In addition this last 
contribution is significant also for equilateral non-Gaussianities, whose effects on the galaxy power spectrum have 
been shown to be negligible, [37] . 

Specifically, for local non-Gaussianities, these two contributions present a large-scale behavior characterized by extra 
factors ofl/k^ and 1/fc^, respectively, when compared to the predictions for the galaxy bispectrum for Gaussian initial 
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conditions. This implies that, for significant quadratic non-Hnearities in the galaxy bias, they are expected to represent 
the leading contribution to the galaxy bispectrum on slightly different ranges of scales, with the component induced 
by the initial trispectrum dominant at the largest scales and the one induced by the initial bispectrum dominant 
on intermediate scales of the order of fc ~ 0.005 /iMpc~^, while Gaussian terms represent the leading component on 
smaller scales. 

For equilateral non-Gaussianities we assumed, for simplicity, a vanishing initial trispectrum. However, the contri- 
bution to the matter trispectrum due to the initial bispectrum alone, still leads to a significant effect on the galaxy 
bispectrum. This appears to be a generic effect of primordial non-Gaussianity, regardless of the model, while no 
analogous effect is present for the galaxy power spectrum. 

We can therefore expect the galaxy bispectrum to be a probe of primordial non-Gaussianity sensitive to the initial 
bispectrum as well as to the initial trispectrum. This is not surprising given the effects that the initial bispectrum has 
on the galaxy power spectrum. However, it is particularly interesting because, despite the fact that the simple local 
model of Eq. ([ij predicts a trivial relation between the amplitudes of the primordial bispectrum and trispectrum, 
with the latter given just by f^i^, in general we can expect the determinations of these two quantities to provide 
distinct constraints on the parameters of the inflationary model. Moreover, the sensitivity of the galaxy bispectrum 
on primordial non-Gaussianity does not seems to be restricted to the local model, as it appears to be the case for the 
galaxy power spectrum. 

A comparison of these results for the galaxy, and halo, bispectrum with numerical simulations is, at this point, 
necessary, to obtain reliable quantitative predictions to be used in the analysis of large-scale observations. Further- 
more, the study of the effect of non-Gaussian initial conditions on the halo bispectrum can significantly improve 
our theoretical understanding of the effect on the halo power spectrum itself, as we can well expect it to point-out 
unavoidable shortcomings and inaccuracies of the different approaches proposed so far in the literature, all assuming 
somehow different approximations. 

Scoccimarro et al. [20^ and Sefusatti and Komatsu |21) have shown that constraints on non-Gaussianity of both 
local and equilateral kind from measurements of the galaxy bispectrum in future redshift surveys are comparable and 
even better than CMB constraints. However, these works assumed the tree-level prediction in perturbation theory 
for the galaxy bispectrum and neglected contributions depending on the matter trispectrum. An update of these 
results is clearly necessary, in light of the results of this paper as well as those of Jeong and Komatsu [3T] and we will 
leave it for future work. If the qualitative results presented here were to be confirmed in numerical simulations, it 
is reasonable to assume that such predictions will significantly improve. In particular, a combined analysis of power 
spectrum and bispectrum of biased populations at large-scales, when an accurate description of non-linearities in the 
bias relation is achieved, might represent the best probe of a departure from Gaussian initial conditions, the latter 
defined by a non-vanishing primordial bispectrum as well as a non-vanishing primordial trispectrum and possibly 
higher-order correlators. 
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